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ABSTRACT 



According to the discussion on Birkhoff's theorem by Peebles (1993 ), a void with a negative per- 
turbation 5 may evolve as a separate homogeneous universe with a local expansion parameter Hy — 
Ho{l — (5/3). This slightly low density "universe" will diverge from the mean, producing a void of ever 
lower density. As a result, the contents of voids will "fall" outward at a velocity equal to the difference 
between the local, and the mean expansion parameters times the radius of the void. Observational con- 
straints on the outfall velocity can be placed as a result of the fortuitous event that void, and non-void 
hydrogen Lya absorbers have distinct characteristics - both in their equivalent width distributions and 
their Doppler parameter distributions - which are clearly distinguished when cloud environ ments are 
measured in terms of the sum of the tidal fields of surrounding galaxies acting on the cloud ( Manning, 



2002| ). These constraints dictate that the radial outfall velocity from an average void of radius ^ 15 h- 



75 

Mpc must be less than 100 kms . The implications of this are probed with a "negative" top- hat 
simulation, using a 1-D Lagrangian code, to determine the relationship between the mass deficit in voids 
and the local expansion parameter in a flat lambda cosmology. The outfall velocity constraint shows 
that the void matter density must be greater than 75% of the mean matter density. I argue that this 
implies that Qm ^ 0.86. Thus the total amount of matter in voids is of order twice the total mass in the 
filamentary structures. 

Subject headings: intergalactic medium — quasars:absorption lines - dark matter 



1. INTRODUCTION 

Voids are structures noted for their pronounced lack of 
galaxies; concentrations of matter are thought to take the 
form of a luminous "filamentary" structure, so that to- 
gether they resemble a "foam" of roughly spherical voids, 
edged by filaments and walls, and sprinkled with richer 
concentrations of galaxies where the filaments intersect. 
While galaxies provide an important tracer of mass in 
the universe, it is possible that there may be significant 
amounts of undetected matter in the form of highly ion- 
ized clouds, or perhaps a significant background of dark 
and gaseous matter, either within the filamentary struc- 
ture, or distributed within the voids. Observations of the 
low-redshift Lya cloud popi ilation have shown that clouds 
cluste r about galaxies {e.g., Tripp et al., 1998 ; Chen et al.. 
2001), but there are also many clouds which appear to be 
isolated. 



In Manning (2002) (hereaf ter Paper 1), a primary cloud 
catalog ( Penton ct al., 2000 ) was separated into comple- 
mentary void and non-void catalogs through the use of a 
parameter which assesses cloud environments by weight- 
ing the effects of galaxy mass and distance. This parame- 
ter, the scalar tidal field, which was summed over galaxies 
within 7.5 h75 Mpc of a cloud, is optimal for this purpose, 
as it may easily be calculated from galaxy catalogs. In ad- 
dition, there is a physical basis for supposing tides affect 
cloud distributions, since it is the tide which determines 
the lower limit of cloud density which is stable to dynam- 
ical disruption. In Paper 1, a dimensionless form of the 
tidal field T was evaluated at the lo cations of each of a 
catalog of low-redshift H I absorbers ( Penton et al. , 20do| ) , 



and used to characterize the clouds as members of a void, 
or non-void environment, according to whether that field 
was greater or less than some limiting tide Tlim- Equiva- 



lent width distribution functions (EWDF) are calculated 
for each sub-catalog (see Eq. 44, Paper 1). These EWDFs 
are well-approximated as linear in \ogdM / dz, vs. logW, 
where TV — n{> Nhi), and W is the rest equivalent width 
(EW) in mA. One may vary Tum, calculate the EWDFs 
of each sub-catalog, and derive the weighted linear fits 
thereof, based on the number of clouds contributing to the 
calculation of the EWDF at a given EW. Ditsributions are 
fitted to the equation, logdN/dz = C 5 log (W/63 mA). 
It was found that the trends of fitting parameters for void 
and non-void catalogs as a function of Tum appear linear 
at large, or small Tum (see Fig. 1), but undergo a non- 
linear change within a relatively small range of T. As it 
happens, there is a corresponding change in the distribu- 
tions of Doppler parameters at a similar range of T for 
clouds sorted by T (see the differential histogram. Fig. 8 
in Paper 1). It appears certain that two distinct popula- 
tions have been isolated, and that a transition zone lies 
within the range, —1.3 <^ logT <^ —0.7. 

Apart from the interest in the detailed nature of these 
clouds, and what they say about the nature of voids (paper 
in preparation) , the very fact that they can be sepated says 
somethin g about the de nsitv of voids. For, by Birkhoff's 
theorem ( |Birkhoff, 19"23|) , if voids have a low density, they 
should behave like a low-density universe, and suffer less 
deceleration, and have larger local expansion parameters, 
than the mean. With simulations I calculate the void ex- 
pansion parameter as a function of the void matter density 
relative to the average. If the smallness of the range of T 
over which the transition of void to non-void clouds occurs 
can be used to constrain the outfall velocity, then this can 
be used to constrain the the void matter density. This is, 
in short, my plan. 

The cosmology assumed for this paper is a standard flat 
lambda model with h — 0.75. The total matter density of 
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Fig. 1. — The trend in slopes (panel a) and intercepts (panel V) 
of equivalent width distribution functions defined by catalogs with 
tidal field lower (top of each panel), or upper limits (bottom), re- 
spectively, for the low-redshift cloud sample (see §6.2.1 Paper 1). 
There is an apparent strong transition in the slopes in the range 
— 1.3 < logT < —0.7 (two dotted virtical lines), which is also seen 
in the intercepts. 



voids is either referred to as fiy = pv / Pcrit, or as /o = 
r2y/r2,„, where — 0.3 is initially assumed. As noted 
in Paper 1, the large discovered line density of void clouds 
appears to require that clouds are to a significant degree 
self-gravitating and discrete. I treat them as such herein. 

2. CONSTRAINTS ON CLOUD OUTFALL VELOCITY 

Peculiar velocities in galaxy clusters cause a phenomenon 
known as the "finger of God" . A similar error in attributed 
position may result when void clouds have large outfall 
velocities. Large outfall velocities will smear the distinc- 
tion between void and non-void clouds, as attributed po- 
sitions are based on the redshift and the mean Hubble 
constant (see §1). Figure 1 shows the trend of linear fit- 
ting parameters for EWDFs defined by various Tum- It 
appears to show a fairly small amount of smearing, sug- 
gesting that the transition zone (dotted vertical lines) is 
—0.13 ^ logT <^ —0.7. How can this observation be used 
to constrain the outfall velocity? The void filling factor as 
a function of provides the means. Tidal fields were 
calculated along the pathlength containing the absorption 
systems (see §5 of Paper 1). Figure 2 shows the fractional 
pathlength with T < Tum- As a measure of the proba- 
bility of a random spot in the universe having a tide less 
than Tum , this is interpreted as the fractional volume with 

< %im- Figures 4 in Paper 1 show that tidal fields in- 
crease rapidly at void edges. By raising the tidal field up- 
per limits, we are increasing the volume of what we think 
of as voids. The bracketing filling factors which correspond 
to the above range of Tum are, 0.75 < /y < 0.91 (see Fig. 
2), with a mean (/y) = 0.86, corresponding to T = 0.1. 

Now, consider an average, fiducial void, which, at Tum = 
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Fig. 2. — The fraction of redshift space with T < plotted 
against logTum- This fraction can be considered to be a fair approx- 
imation of the volume filling factor for voids defined by T < Tjim- 
Dotted lines refer to the region of transition similarly highlighted in 
Fig. 1. 



0.1, has a radius of 15 Mpc (~ 11.25 h"! Mpc) (e.g., 
Lindner et al., 1995). The range of uncertainty in the 



position of the transition zone corresponds to radii of 
rv = 15(/y/0.86)i/3 = 14.33, and 15.396 Mpc, for void 
filling factors fv = 0.75 and 0.91, respectively; an average 
variation of ±530 kpc. This uncertainty in the true fidu- 
cial radius corresponds to a line of sight (LOS) velocity 
error of Hq 5r ±40 kms^^, and would roughly repro- 
duce the same range of uncertainty apparent in the transi- 
tion zone. Admittedly, many clouds will be moving along 
vectors at large angles to the LOS, and so even with high 
outfall velocities, they could have low LOS peculiar veloc- 
ities. But if clouds are randomly distributed in the void 
relative to our LOS, and radially outfalling, an average of 
half of the clouds would have to be moving at an angle less 
than 60 degrees from the LOS, so that the observed LOS 
outfall velocity at the edge would be greater than half of 
the radial outfall velocity for half of these clouds. However, 
if clouds and mass are uniformly distributed in voids, then 
half of the clouds will be at a radius r > 0.794 ry, so that 
most clouds would have an outfall velocity Vout > 80% 
of the maximum. Therefore, we expect that half of the 
clouds will have a peculiar outfall velocity along the LOS 
^ 0.4 times the radial outfall velocity at the void edge. 
Since the observed range of the LOS positional error cor- 
responds to a velocity error of ~ 40 kms^^, this would 
seem to limit the true radial outfall velocity at the void 
edge to be Vout ^ 40/0.4 = 100 kms^^. However, since we 
know that some of the error in the determination of the 
transition is plausibly due to peculiar velocities of non-void 
clouds in proximity to galaxies (see §5.2, Paper 1), this is 
probably an over-estimation of the upper limit on Vout- 

3. THE OUTFALL VELOCITY FROM VOIDS 

Voids are thought to grow more rapidly th an the scale 
factor {e.g., Bcrtschingcr, 1985 ; Piran, 1997 ). This can 
be conceptually understood by comparing the evolution of 
the expansion parameter for a flat (the overall universe), 
and an open universe (the void). For the former. 
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and the latter, ( ^cott et al., 2000| ), 



Hv{z) - Hv{0)Vifo ^mz + l- f^A)(l + zy + n^, (2) 

where Hy reflects a Hubble constant in the void, and /o 
is defined in §1. One might think that as z approaches 
infinity, their respective expansion parameters would ap- 
proach the same value so that the ratio of their current 
values would simply be. 



nrnii + z)^ + nA 



{fon^z + i-nA)ii + zy + nA 



1/2 



(3) 



so that a "peculiar" outfall velocity could be determined 
by subtracting off the mean Hubble flow. 



(7^H-l)i^ory, 



(4) 



where ry is the radius of the void. The dashed line in Fig. 
3 shows the functional form of Eq. 4 as a function of /q. 

However, while the theorem says that the lower-density 
void acts like a separate "universe", our reconing of it, in 
terms of the correspondence between lookback time and 
redshift, is not straightforward, for the lookback time to 
a given redshift is a function of cosmology, and a flat uni- 
verse is younger than an open universe with the same Hq: 
we cannot so easily derive our answer. 

However, we may confidently approach this subject by 
simulating the void as a "negative" top-hat with pertur- 
bation S. A 1-dimensional Lagrangian hydro/gravity code 
(Thoul & Weinberg, 1995) was acquired courtesy of one of 
its authors, (A. A. T.), and adjusted for a fiat, A cosmol- 
ogy (see §1). A more detailed description of the code and 
some of its broader uses will be presented in a subsequent 
paper (Manning, in preparation). Linearized analysis of 
voids dictate that the appropriate initi al velocity perturba- 
tion should be v 

Peebles, 1993), where Hy is the local 
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S/3)Hv{z)r (Bertschinger 

void expansion pa- 
rameter. Simulations of the dark matter component were 
performed using this velocity perturbation, with a range 
of uniform underdensities S embedded in the otherwise fiat 
universe, beginning at z = 50. A 5 = case was run to 
normalize Hq. The expansion parameters at z = were 
derived. The outfall velocity was calculated by subtract- 
ing Hq from Hv, and multiplying the result by the 15 
Mpc of our fiducial void. Figure 3 shows the results (solid 
line) , where Vout is the outfall velocity. For void densities 
/o <, 0.2, outfall velocities are Vout > 370 kms~^. How- 
ever, this analysis (§2) finds that Vout ^ 100 kms^^, and 
thus, /o ^ 0.75. We may apparently conclude that the 
Hubble fiow in voids is very close to the average Hubble 
fiow, and that therefore the difference in density between 
voids and the mean is minimal. 

4. IMPLICATIONS FOR Qm 

Using the mean filling factors fv — 0.86 for voids and 
fp = 0.14 for filaments, we may analyze the distribution 
of matter implied by the above. The basic equation for 
this analysis is. 



^771 — fv fo + If ^f, 



(5) 




Fig. 3. — The peculiar outfall velocity as a function of void density 
relative to the mean for a void of radius ry = 15 Mpc, and h = 0.75. 
Results for a range of underdensities S are shown as filled circles for 
f^m = 0.3, and as open circles for Qrn = 0.85 (see §4). The dotted 
lines refer to characteristic values of fo noted in the text. The 
dashed line represents the outfall velocity as derived from Eqs. 3 
and 4, which ignores the disparity between time-scales in flat and 
open cosmologies. 



where faflm is substituted for Qy- 
ftm — 0.3 and fo > 0.75, we find. 



Solving for Qp, when 



^m{l-fvfo) 
fF 



0.761, 



(6) 



which is absurd, since gravitational clustering could not 
occur unless fij? > 1.0. It is apparent that Clp must be 
much greater. Introducing commonly accepted values into 
Eq. 6, for instance, il™ — 0.35 an d fp — 0.18, the mean 



of the range 0.1 ^ fp 0.26 {e.g., |Cen fc Ostrikcr, 1999 
El- Ad fc Piran, 1997]) , we find ^ 2.2. If we assume 



this more reasonable value of ilp, then it turns out that 
using the above constraint on /o, then 



fp Hp 
l-fvfc 



;> 0.86. 



(7) 



That is, a low outfall velocity from voids mandates a 
"high" But is the outfall velocity for large fi^ the 

same function of /o? When the 1-D code was run for 
f2„i = 0.85, a graph consistent with that for = 0-3 re- 
sulted (dashed line, open points. Fig. 3), so the same con- 
straints maintain in a high Qm universe as with flm = 0.3. 

5. DISCUSSION AND SUMMARY 

5.1. Impact on Filaments 

First, let us consider the effects of the deposition of 
energy in volved in the outfall of void clouds onto fila- 
ments. If clouds are diffuse, and the kinetic energy of 
outfall is transformed entirely into random motions, then 
T = /xm//UQ,jj/3fc, where k is the Boltzmann constant, 
and fi ~ 0.59 is the mass of an average particle. A large 
velocity Vout = 400 kms~^ would produce a plasma of 
temperature nearly 4 x 10^ K, similar to that predicted by 
Cen & Ostriker (1999) for filaments, suggesting that ab- 
sorbers should have Doppler parameters b ^ 250 kms^^. 
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much larger than observed 6- values of Lya lines of non-void 
clouds {b w 60 kms~^, Paper 1). For Vout — 100 kms~^, 
diffuse clouds are converted to a temperature T = 2.4 x lO'^ 
K; 6 ~ 62 kms^^. Interesti ngly, Far Ultravio let Spectro- 
scopic Explorer observations ( phuU et al., 2000 ) find actual 
6— values may be roughly half that of the observed Lya 
line, the excess being attributed to bulk motions. Doppler 
parameters of order 30 kms~^ may be explained if clouds 
are not diffuse, but held by dark matter halos, and hence 
centrally condensed. In this case, only gas with a density 
of order that of the filament will be stripped, l eaving the 
centr al cloud to p roceed relatively unmolested ( Murakamj 
& Ikepichi, 1994), minimizing the t hermal broaden ing of 



the cloud. Shear-induced vorticity (Manning, 199£), will 
contribute to the Doppler broadening of the cloud, and 
ram pressure will heat it, so that these absorbers should 
have greater Doppler parameters than void clouds, as ob 
serve d (Pap er 1 
ning [1999| ) 



§5.6). According to the scenario in |Man 
the induced vorticity of clouds, which are 
falling into dissipated gaseous envelopes of galaxies and 
groups of galaxies, introduces an orderly cycle of com- 
pression and cooling which results in elevated central H 
I equivalent widths as the cloud approaches mass concen- 
trations. This inverse correlation between cloud rest equiv- 
alent widt hs and projected galactocentric radii ha s been 
obser ved (Lanzetta et al., 1995 ; Tripp et al., 1998; Chen| 



et al., 1999; Chen et al., 2001) 



5.2. Implications for Cosmology 

The conclusion that voids have a high matter density 
depends heavily on the successful separation of void, from 
non-void clouds. That the Doppler parameter differential 
histograms, and trends in EWDF fitting parameters both 
implicate the same transition zone suggests that this is a 
robust result. Outside the transition zone — 1.3 < logT ^ 
—0.7, the trends seen in Fig. 1 suggest distinct cloud 
characteristics with few interlopers. Within this range, 
however, the trends change in a manner consistent with 
the effects of peculiar velocities of non-void clouds about 
galaxies, compounded with that of the outfall of clouds 
from voids with Vout <, 100 kms^^ {Hy <^ 1.1 Hq). The 
simulations have shown that for Vout <, 100 kms^^, the 
ratio of the void density to the mean is /o ^ 0.75, requir- 
ing to be of order 1. Under the strictures of the flat 
Lambda model, this requires SIa « 0, and A « 0. 

The current standard model postulates two forms of pos- 
itive energy which, under the analysis of the multipole fluc- 
tuati ons in the CMB, appear t o imply that + — 1 
{e.g., ie Bernardis et al., 200C). Independent assessments 
of flm are all based on concentrations of matter, usu- 
ally involving dynamical tests, or biasing estimates {e.g. 
Fukugita et al., 1998; Verde, 2001), and result in estimates 



rim ~ 0.3. However, since these methods are insensitive 
to a possible background, or additional quantities of mat- 
ter in void regions, these estimates are necessarily under- 
estimates. As we have seen, constraints on the velocity of 
outfall from voids imply a rather large "background" . 

On the other hand, supernova la studies str ongly imply 
an exponential expansio n of the universe {e.g., Ricss et al.J 
1998; Filippenko, 2001). How is this to be reconciled? 
What has not been shown by observation or theoretical 
derivations is that what is created in conjunction with the 



exponential expansion, and what enables the universe to 
be flat, is a "dark" energy. The conclusion that ilm ^ 0.85 
and A ^ can only mean that a dark energy is not formed, 
but in its stead, ordinary matter, dark and baryonic. 

It is unclear how one might adjust the parameters of 
the void evolution simulations to be consistent with these 
conditions, for if matter is forming in voids, the tendency 
of the mass density to rapidly decline is reduced. A careful 
numerical modeling of void cloud absorbers as remnants 
of sub-galactic perturbations (paper in preparation) may 
shed light on this problem (Manning, in preparation). 
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#AST-0097163 and the UC Department of Astronomy. 
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